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We present explicit closed-form expressions for the two-loop Euler-Heisenberg 
Lagrangians in a constant self-dual field, for both spinor and scalar QED. The 
simplicity of these representations allows us to examine in detail the asymptotic 
properties of these Lagrangians, and to construct their imaginary part using Borel 
dispersion relations. In particular, for this self-dual case we obtain the explicit 
form of the Lebedev-Ritus functions appearing in the Schwinger representation 
of the imaginary part at two loops. Using the connection between self-duality 
and helicity, we also obtain explicit formulas for the low energy limits of the 'all 
+ helicity' N photon amplitudes, in scalar and spinor QED at one and two loops. 

1 Introduction 

In gauge theory a very special role is played by self-dual fields, i.e. fields 
satisfying the condition 



Fields of this type are prominent in QCD for a number of different reasons: 

• Instantons are self-dual. 

• Among all covariantly constant gluon backgrounds, only the self-dual 
quasi-abelian background is stable (at one-loop) under fluctuations [p]]. 

• Large classes of integrable models can be obtained by dimensional re- 
duction starting from self-dual Yang-Mills theories 0. 

• Self-dual fields are helicity eigenstates, so that the effective action in such 
a field carries the information on the corresponding gluon amplitudes 
with all equal helicities 0. Such 'all +' amplitudes generally exhibit a 
particularly simple structure, at the tree level || and beyond []5J. 
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In the abelian case most of this motivation does not exist, with the exception of 
the last point mentioned, which holds in the abelian well. Correspond- 

ingly, little use has been made so far of self-dual fields in QED. In the present 
contribution, we will consider the case of constant self-dual fields in QED, and 
show that the self-duality condition leads to very remarkable simplifications 
for the effective action in such a background field ||. 



2 One-loop Euler-Heisenberg Lagrangians 

Before presenting these two-loop results, let us shortly recapitulate some facts 
about Euler-Heisenberg Lagrangians. At one-loop, the on-shell renormalized 
effective Lagrangians in a constant background field, for spinor and scalar 
QED, are given by the well-known formulas ]7], |8| 
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Here a, b are related to the two invariants of the Maxwell field by a 2 — b 2 = 
B 2 -E 2 ,ab = E- B. 

These effective Lagrangians are real for a purely magnetic field, while in 
the presence of an electric field there is an imaginary (absorptive) part, indi- 
cating the process of electron-positron (resp. scalar-antiscalar) pair creation 
by the field. For example, in the case of a purely electric field, E, the effective 
Lagrangians (|2.2| ) have imaginary parts given by: 
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where f3 = eE/m 2 . These expressions are clearly non-perturbative in terms of 
the field and coupling. Their physical interpretation is that the coefficient of 
the fc-th exponential can be directly identified with the rate for the coherent 
production of k pairs by the field M. 
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3 Two-loop Euler-Heisenberg Lagrangians 



The two-loop corrections to the Euler-Heisenberg Lagrangians, taking into 
account an additional photon exchange in the loop, where first calculated by 
Ritus in the seventies (see also []nj). More recently these Lagrangians have 
been recalculated [|ll|] using the 'string-inspired' formalism |P2| , 13 , | 



l r 



However, in all cases the results are not nearly as explicit as the one-loop 
formulas (|2.2| ); they involve two-parameter integrals and a counterterm from 
one-loop mass renormalization. As far as the magnetic case is concerned, where 
the effective Lagrangian is real, the complicated form of these representations 
is perhaps not too bothersome. For values of the magnetic field small compared 
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to the 'critical' field strength B c = — , the Lagrangian can be computed using 
the weak field expansion, whose coefficients are easy to compute to fairly high 
orders 
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For general values of the field strength numerical integration can be used. 

Calculating the imaginary part of the corresponding Lagrangian for the 
electric field case is a more difficult matter. At the one-loop level, the imag- 
inary parts ( |2.3| ) can be obtained from ( |2.2j ) by a simple application of the 
residue theorem. The analogous analysis for the two-loop parameter integrals 
is already highly nontrivial. Nevertheless, Ritus and Lebedev [T7| were 
able to obtain along these lines the following two-loop generalization of the 
Schwinger decompositions (|2.3j ): 
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where a = j- is the fine-structure constant. The coefficient functions _K"^ pm,scal (^ 
appearing here were not obtained explicitly by fl7 |. However, it was shown 
that they have small (3 expansions of the following form: 
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Note that to the order given K^ m = i^| cal . For k > 2, these expansions start 
with terms that are singular in the limit of vanishing field (3 — > 0, which seems 
to be at variance with the fact that these coefficients have a direct physical 
meaning. In [17 a physically intuitive solution was offered to this dilemma. 
Its basic assumption is that, if one would take into account contributions from 
higher loop orders to the prefactor of the fc-th exponential, then one would 
find the two lowest order terms in the small - (3 expansion of l^ pm ' scal (/3) to 
exponentiate in the following way, 
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It would thus be possible to absorb their effect completely into a field-dependent 
shift of the electron mass: 



m*(k,E) = m + -akckVeE — -akeE /m 



(3. 



Moreover, these contributions to the mass shift have a simple meaning in 
the coherent tunneling picture [TTj: The negative term can be interpreted as 
the total Coulomb energy of attraction between opposite charges in a coherent 
group; the positive term, which is present only in the case k > 2, represents the 
energy of repulsion between like charges. For the attractive term, a completely 
different derivation of the same exponentiation was given in [18 



4 The imaginary part via Borel dispersion relations 

Clearly, it would be of interest to understand in greater detail this prefactor 
series i^ pm ' scal (/?). In an effort to learn more about it, in [|TjJ we used Borel 
techniques to study the imaginary part of the two-loop effective Lagrangian 
for a constant electric field background. 
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Let us review some basic facts on Borel summation. Consider an asymp- 
totic series expansion of some function f(g) 

oo 

f(g)~J2*n9 n (4.9) 

n=0 

where g — > + is a small dimensionless perturbation expansion parameter. In 
many physics applications perturbation theory leads to a divergent series in 
which the expansion coefficients a n have leading large-order behaviour 



a n ~ (-l)VT(/in + z/) (n->oo) (4.10) 

for some real constants p, p > 0, and v. When p > 0, the perturbative expan- 
sion coefficients a n alternate in sign and their magnitude grows factorially, just 
as in the Euler-Heisenberg case. Borel summation is a useful approach to this 
case of a divergent, but alternating series. The leading Borel approximation is 




(4.11) 



For a non-alternating series, we need f(—g). The Borel integral Q4.11D is 
an analytic function of g in the cut g plane: |arg(g)| < n. So a dispersion 
relation (using the discontinuity across the cut along the negative g axis) can 
be used to define the imaginary part of f(g) for negative values of the expansion 
parameter: 
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Returning to the Euler-Heisenberg Lagrangian, for a uniform magnetic 
background the weak-field expansions of the one-loop Lagrangians ( |2.2[ ) are 
precisely of the form ( f4.9|) , (|47TO|) with g = (^t) 2 . For example, in the spinor 
QED case one has 
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where the B m are Bernoulli numbers. For a uniform electric background, the 



only difference perturbatively is that B 2 is replaced by — E 2 ; that is, g = (f^l) 2 
is replaced by — g = — (f^) 2 - So the perturbative one-loop Euler-Heisenberg 
series becomes non-alternating. Then from ( |4.12|) , with p = \ and /j, — v — 2, 
we immediately deduce the leading behaviour of the imaginary part of the one- 
loop Euler-Heisenberg effective Lagrangian: 
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Taking into account also the sub-leading corrections ( |4.13| ) to the leading large- 
order behaviour of the expansion coefficients a„ , one can apply ( 4.12 ) succes- 
sively to reconstruct the full Schwinger series (|2.3| ). 

At the two-loop level, no closed-form expression is known for the corre- 

(2) 

sponding expansion coefficients a n , so that the corresponding analysis be- 
comes much more involved. In [IB], first the integral representation for the 
purely magnetic case obtained in [|TTj] was used to compute these coefficients 
up to n = 15. It was then established by a numerical analysis that 
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Thus the leading large-order growth corresponds to the form ( (4.10p , and more- 
over differs from the one-loop case (|4.13| ) only by a global prefactor. To the 
contrary, the subleading correction term given in ( (4.15[ ) shows a much weaker 
n dependence than is found for the first correction in the one- loop case ( |4 . 1 3| ) . 
This means that in the two-loop case the dominant corrections are to the 
prefactor in the leading behaviour. This is in contrast to the one-loop case 
where the first correction to the leading behaviour is exponentially suppressed. 
Indeed, applying the Borel relations, the correction term ( |4.15| ) leads to 
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Thus the structure of Q4.16 ) conforms already to the form ( |3.5|) ,(| 
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5 The self-dual case 



It seems very difficult to make further progress along these lines for the purely 
magnetic/electric field cases, let alone the general constant field case. To the 
contrary, in the self-dual case this analysis can be carried much further [HJ. 
Technically, this is because the self-duality condition ( |1 . 1|) implies that 



F = -f 2 l (5.17) 

where f 2 = \F ia ,F tlv , and 1 denotes the identity matrix in Lorentz space. 
This leads to enormous simplifications for this type of calculations. 

In the self-dual ('SD') case, the integral representation (|2.2| ) for the renor- 
malized one-loop scalar effective Lagrangian becomes 



scal { ' (4vr) 2 4*; 2 J Q t 3 

2 

where k = - ^= is the natural dimensionless parameter. The Lagrangian for 

the spinor QED case differs from this only by the standard global factor of —2 
for statistics and degrees of freedom: 
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4r D V) = -2£«rV) (5.19) 

Note that this relation holds for the renormalized effective Lagrangians, not 
for the unrenormalized ones. It is due to a supersymmetry of the self-dual 
background [111. For real k, this self-dual Lagrangian is real and has properties 



very similar to the magnetic Lagrangian. Similarly, for purely imaginary k the 
self-dual Lagrangian has an imaginary part, and provides a good analogue 
of the electric field case. These cases will therefore be called 'magnetic' and 
'electric' in the following. In particular, the imaginary part of the 'electric' 
Lagrangian has a Schwinger-type expansion 
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Surprisingly, for the self-dual case even at two loops all parameter integrals 
can be done in closed form, leading to the following explicit formulas [0]: 
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Here we have introduced the function £, 



f (x) = -x(i){x) - ln(x) + — ) (5.22) 

and i/j(k) = i^l n r( K ) i s the digamma function. This function has very simple 
expansions at zero as well as at infinity, 



1>(x) ~ ---T + Et- 1 )^)^ -1 (5-23) 

' k=2 
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so that the formulas ( |5.21| ) directly yield closed-form expressions for the co- 
efficients of both the weak-field and strong-field expansions of the self-dual 
Lagrangians. Moreover, it is easy to obtain also a closed representation for 
their imaginary parts which exist in the 'electric' case: 
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Note that these formulas display the same structure found by Ritus and Lebe- 
dev for the physical electric case, (|3.5|),(|3~^), but that in the self-dual case we 
have simple closed formulas for all the prefactors of the Schwinger exponentials. 
In ( |5.25p we have written these prefactors as series in -, which corresponds 
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to the (3 - expansion (|3.6j). It is easy to see that this expansion is asymp- 
totic, rather than convergent, for all k. For the scalar QED case, we have 
also verified that the first (k = 1) exponential in ( |5.25| ), including its com- 
plete prefactor, can be correctly reproduced by an analysis of the asymptotic 
behaviour of the weak field expansion coefficients and the application of the 
Borel dispersion relations. It is also immediately evident that, for any given 
value of ft, the two-loop contribution to the prefactor of the k-th exponential 
will dominate over the corresponding one-loop quantity in ( |5.20|) if k is taken 
large enough. Given the close similarity between the electric and 'electric' 
cases at the one-loop level it is natural to assume that these properties hold 
also true for the electric case. 



6 The N - photon amplitudes 

As is well-known, the Euler-Heisenberg Lagrangians ( p.2|) can be used to obtain 
the one-loop QED N - photon amplitudes in the low-energy approximation, 
i.e. for photon momenta fcj such that m 2 is much larger than all ki ■ kj 's. 
After specializing to a self-dual background it still contains the information 
on the component of the iV - photon amplitude with all helicities equal, say, 
all '+'. Using (|2.2|),( ]oTi2"T| ),( |5.24| ) together with the standard spinor helicity 



formalism (see, e.g., |2(J), one obtains for this low-energy limit the following 
explicit formulas: At one loop, 
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and at two loops 
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In these formulas, the information on the external momenta is all contained 
in the invariant Xn, 



XN = i^lJ [i2] 2 [34] 2 ---[(iV- 1)A^] 2 + all permutations [ 
2 2 ^ J 

Here [ij] = (kf\kj) denote the basic spinor products ]20 



7 Conclusions 

The existence of the simple closed-form expressions ( |5.21| ) for the self-dual two- 
loop Euler-Heisenberg Lagrangians is surprising, and we do not know of any 
comparable result in gauge theory. These results have allowed us to perform 
a more complete analysis of these effective Lagrangians than is possible for 
other backgrounds |J, including the complete elucidation of the structure of 
the imaginary parts. Moreover, we have used them to obtain explicit formulas 
for the low-energy limits of the 'all +' components of the two- loop iV - photon 
amplitudes. Let us close with a remark on the behaviour of these amplitudes 
in the limit where the number of photons becomes large. It is easy to show 
that 
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The Borel dispersion relations allow one to relate this quantity to the k — ► oo 
limit of the corresponding ratio of the imaginary parts ( |5.20[ ),( |5.25| ), and thus 
to the factor air appearing in the Lebedev-Ritus exponentiation formula ( |3.7|) . 
Based on the expectation that this exponentiation occurs also in the self-dual 
case, we conjecture that at loop order I, the weak-field expansion will take the 
form 

„ K) J0_^ C ^ ( , 29) 
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where for each I, the expansion coefficients cl have the same leading large n 
growth rate This would imply that 
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It would be very interesting to verify this identity at the three-loop level. 
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